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Introduction 



The radiative kaon decay amplitudes K + — > n + e + e~ (fi + are of great interest of the chiral 
perturbation theory (ChPT) [HE} [3] because in the lowest order of ChPT the decay amplitudes 
are equal to zero pQ EJ El [7J, E] . There two opinions about the next order of ChPT. 

The first opinion is that in the next order of ChPT the baryon (or quark) loops dominate [H, 
E]. As these fermion loops also determine meson form factors in the low energy region [T], in this 
case, one can supposes that ChPT calculation [H El El H] points out possible relations between 
the low energy parameters of meson form factors and the differential rates of the radiative 
kaon decays. These relations can arise if we keep in the ChPT diagrams the real vector meson 
propagators and take into account the quantum numbers of the nearest resonances in possible 
vertices pfTO]. 

The second opinion is that in ChPT the meson loops dominate. It was shown [HI [71 E], 
in the framework of the accepted approach to Standard Model [HI H21 H3] with the point-like 
approximation of weak interactions, that these meson loops can completely destruct the meson 
form factor structure of the radiative kaon decay amplitudes in the low energy region. 

In this paper, we show that the situation with ChPT for radiative kaon decay amplitudes 
is more complicate. There are two types of the meson loops. The first of them are provided by 
the normal ordering of the weak static interactions, and the second are retarded ones. 

Recall that static interactions arise in the Hamiltonian approach to the Standard Model 
(SM) of electroweak (EW) interactions [HI [15] in contrast to the conventional one [HI [[2] based 
on heuristic Lorentz gauge formulation [HUES], where the static interactions are absent. These 
static interactions suppress any retarded meson loop contributions [3, [7J E] that can destruct 
the meson form factor structure of the decay amplitudes. The meson loop contributions can 
be only the tadpole loop diagrams following from the normal ordering of the static interaction. 
This ordering results in an effective action with AT = 1/2 rule [UJ E5] with one unknown 
parameter g$ that can be fixed from other decays as g$ — 5.1. The dominance of weak static 
interactions justifies the application of low energy chiral perturbation theory [U EJ E] as an 
efficient method of description of kaon decay processes [H El l6l [T8]. 

In this paper we study the possibility of extracting information about the meson form factors 
from the K + — > it + e(n) + e(n)~ processes at the current level of the experimental precision. 

The plan of the paper is the following. In Section 1 we present the explicit expressions 
for amplitudes of the processes K + — > in terms of meson form factors. In Section 2 

we discuss possibilities of the corresponding experimental tests. Manifestations of the static 
interactions in decay rates in SM are discussed in Section 3. 

1 Relations between form factors and radiative K decay 
amplitude in ChPT 

1.1 Chiral bosonization of EW interaction 

It is conventional to describe weak decays in the framework of electroweak (EW) theory at the 
quark QCD level including current vector boson weak interactions [HI H2] 

where = £^7^(1 —75)^, d! = d cos 9c + s sin 9c, and 9c is the Cabbibo angle (sin 9c = 0.223). 
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However, a consistent theory of QCD at large distances has not been constructed yet. 
Therefore, the most efficient method of analysis in kaon decay physics [U El El HH] is the 
ChPT [2, 3J. The quark content of vr+ and K + mesons tt+ = (d,u),K + = (s,u),K° = (s,d) 
leads to the effective chiral hadron currents 4 in the Lagrangian (pQ) 



4 = [4±*4] cosfl c + [4±*4] sm8 G , 
where using the Gell-Mann matrices X k one can define the meson current as [2] 

i E A "4 = iXk (K - Kf = Fy%e-*, 
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In the first orders in mesons one can write 



V~ = V2 [smd c (K~dy - Ti'd^K-) + C os6 c (t^tt - n°d^-) } + ... (5) 

and 

A- = V2F n {d?K- sin^c + cos6 c ) + (6) 
here F n ~ 92.4 MeV. The right form of the chiral Lagrangian of the electromagnetic interaction 
of mesons can be constructed by the covariant derivative d^x^ ~ > D^x^ = ^ ^AOx*) 
where = K \ 7r± - 

We suppose also that the quark content of the mesons determines hadronization of QCD 
[T9l [20] conserving its chiral and gauge symmetries. 



1.2 The K+ -> TT+l+l- amplitude 

The result of calculation of the amplitude of the process K + — > (I = e; /i) in the 

framework of the chiral Lagrangian (Q])-(E1) including phenomenological meson form factors 
denoted by fat dots in Fig. Q] takes the form 



T, 
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2g s eG EW L u D^ ( g ) + Pfi ) T (q 2 , k 2 ,p 2 



where g% = 5.1 is the effective enhancement coefficient [5l \T7\. 
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is the coupling constant, L M = IjJ is leptonic current, 
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and f^^(q 2 ) are meson form factors. On the mass-shell the sum ((9l) takes the form 
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Figure 1: K + — ► 7r + 7* diagrams for effective Lagrangian. 



The K + — > 7T + / + / _ amplitude vanishes at tree level [U [5], where the form factors are equal to 
unity: T(q 2 )\jv = jA =1 = 0. 

In terms of the two standard Dalitz plot variables q 2 and t 2 representing the squares of 
invariant masses of l + l~ and pairs, respectively, the amplitude ( fTtJl) leads to the following 
decay rate for the transition K + — > 



(M K - m7T ) 2 



t 2 max (l 2 ) 



T(q\t 2 )=C J dq 2 .\F{q 2 )\ 2 j dt 2 • V (q 2 ,P 

4m? 



(11) 



where 
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(2t 2 + q 2 - 2ml ~ 2mf)(2Af|- + 2mf - 2t 2 - q 2 ) + 
+q 2 (q 2 - 2M 2 K - 2ml) 



(12) 



and 



F{q 2 



(4tt) 2 T(^ 



2 V f^q 2 ) + f^ 2 



mz 



Ml -m 2 



(13) 



The ^-dependence does not contain information about the combination of the form factors 
F(q 2 ), which is of our interest. Integration of (fill over t 2 yields 



(M K -m„) 2 



r(q 2 )=C I ^Lp(q 2 )\F(q 2 )\ 2 . 

4mf 



(14) 



Here (see [5]) 
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(15) 



P(g 2 ) = (l-^) V2 (l + ^) A 3 / 2 (l,g 2 /Mi,m^/Ml), 
X(a, b, c) = a 2 + 5 2 + c 2 — 2(a6 + be + ca), 
and S1C1C3 is the product of Cabibbo-Kobayashi-Maskawa matrix elements V u dV us . 

2 Parameterization of F{q 2 ) 

2.1 Parameterization with meson loops 

It follows from (1T3~|) . taking into account f^{q 2 ) ^ fxiq 2 ) and f^{q 2 ) ^ fxi^ 2 )'- 

F{q 2 ) = [/ V) - /V)] • (16) 

We discuss the differential — > decay rate ( 1TT1 fT4l ) in the ChPT [U [2] with pion and 

baryon loop contributions leading to meson form factors [US]: 



f v { q 2 ) = i + M p -V + « n.(g 2 ) + ... , 

fV) = 1 + M-V + - 1 J 



We parameterize the terms linear in q 2 (determined by the baryon and meson loops [Tl[4"l[6l[T8]) 
by the values of resonance masses [22] M p = 775.8 MeV, I G (J PC ) = 1 + (1~) and M a = 984.7 
MeV, I G {J PC ) = l-(0++), 

4 77T, 2 

°» = 3 ■ wk? = 01926 ' (18) 

and the nonlinear term of the pion loop contribution [T], [2j is given by 

1/2 / W 2 \ - t 



n.(t) 



, ^ /1 \ i/2 / t 1 ' 2 \ - t 

^(jV arctan l(I^J -1 ' * = (W <1; 



t-i/ i\ 1/2 r , ?/2 + (t_ 1)1/2 



(19) 



In order to introduce the resonant behavior of the form factors, the following Pade-type ap- 
proximations [21] to the expressions (fTTj) are considered: 



fftf) = 7[i-{M-V + «d-n,r(g a )}/7]" 1 + (i-7), ram 

AV) = (1 - M-V^" 1 



Here the parameter 7 = 0.85 effectively accounts for higher order loops, and is chosen such as 
to put the position of maximum of fY(q 2 ) to q 2 = M 2 . 

2.2 Predictions for integrated and differential decay rates 

The form factors (1201) lead to the following decay branching ratios and muon/electron ratio 
R = Br(K+ -> n + ii + fj 1 -)/BT{K + -> vr+e+e-): 

Br(K + -> vr+e+e") = 3.88 x 10" 7 , Br(K + — > it + n + n~) = 1.23 x 10~ 7 , # = 0.318. 
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z=q 2 /M 2 z=q 2 /M 2 



Figure 2: Differential rates of (a) K + — > n + e + e~ and (b) K + — > 7r + / u + / u _ decays as functions 
of 2 = q 2 /M^. Solid lines: F(g 2 ) parameterization (1201) : dotted lines: differential distribution 
measured for K + — > 7r + e + e~, z > 0.1 by [23] and extrapolated for z < 0.1 and for — > 
using a model [8j. Given the large sensitivity of our calculation to values of M p and 
M a , the calculation agrees with the experimental data. 

These branching fractions are highly sensitive to the values of M a and M p used in the param- 
eterization: 

ABr(ee)/Br(ee) « 12(AM /M ), ABr(/^)/Br(/^) « 10(AM o /M o ), , . 

ABr(ee)/Br(ee) w -20(AM p /M p ), ABr(/^)/Br(/^) » -17(AM„/M p ). [ 1 

However this sensitivity largely cancels in the muon/electron ratio: 

AR/R » -2(AM /M ), AR/R « 2(AM P /M P ). 

The sensitivity to the parameter 7 is smaller than to the resonance masses: 

ABr(ee)/Br(ee) » -0.5(A 7 / 7 ), ABr(/x//)/Br(///x) w -0.9(A 7 / 7 ), AR/R « -0.4(A 7 / 7 ). 

Taking into account that the relative uncertainties of resonance masses M a and M p that are 
about 1%, our predictions can be roughly quantified as follows: 

Br(if+ -> 7T + e + e _ ) = (3.9 ± 0.8) x 10~ 7 , 
Br(K + -> 7r~ t > + /0 = (1-2 ±0.3) x 10~ 7 , 
i2 = 0.32 ±0.01. 

Differential rates of K ± — > ir^e + e~ and — > 7r ± /i + /i~ decays corresponding to the pa- 
rameterization ( 1201 ) are presented in Fig. [2] along with the rates calculated extrapolating the 
available experimental data on K ± — > 7i ± e + e~ decay [23] using a model [8]. Note that the ex- 
perimentally accessible kinematic region of the K + — > n + e + e~ decay is limited by a condition 
z = q 2 /M] c > (M^o/Mk) 2 « 0.08, while for the — > 7r ± / u + /i~ decay the whole kinematic 
range is accessible. 

The function \F(q 2 )\ determined by the relations ( 1T61) . (1201 ) is presented in Figj3l Its shape, 
being approximated in terms of a linear form factor F(q 2 ) = Fq(1 + A • q 2 /M] { ), leads to a form 



6 




factor varying from a minimum of A = 1.4 for z = q 2 /M] ( = to a maximum of A = 4.9 for a 
point z = q 2 /M] c = 0.32 corresponding to q 2 ~ M 2 . An effective average form factor slope of 
the K ± — > TT ± e + e' decay as would be measured by an experiment in the accessible kinematic 
region q 2 > M 2 is estimated to be A ~ 2.3. This value should be subject to variation among 
different experiments, depending, in particular, on the experimental acceptance as a function 
of q 2 . 



2.3 Prospects for experimental tests 

The experimental data for branching fractions, their ratio, and the form factor slope (not 
considering a slope measurement in the K ± — > 7r ± y u + / u _ channel which is subject to large 
uncertainties) are [22] : 

By(K + -> 7r + e + e~) = (2.88 ± 0.13) x 10~ 7 , 
Br(K + -> 7T+//+// - ) = (0.81 ± 0.14) x 10" 7 , 
R = 0.281 ±0.050, 
A = 2.14 ±0.20. 

The experimental precision of Br(K + — > n + e + e~) is mainly determined by a single measure- 
ment |23| . The experimental uncertainty of Br(K + — > -K + fi + fi~) is dominated by a PDG error 
scale factor |22| emerging from inconsistency of three measurements [241 1251 [26] . Significant 
experimental improvements are expected in near future, when the data sample collected by the 
NA48/2 experiment at CERN is analyzed. 

Our predictions for branching fractions of the two decays, their ratio R, and the effective 
form factor slope are in agreement with the experimental data. On the contrary, meson domi- 
nance models [27j fail to describe the effective form factor slope, predicting substantially lower 
slope values. 

Monte Carlo simulations involving realistic estimations of experimental conditions show 
that a deviation of K ± — ► 7r ± e + e~ event distribution predicted by (1201) from a distribution 
corresponding to a linear form factor F(q 2 ) (i.e. the predicted dependence of the effective form 
factor slope A on q 2 ) can be experimentally detected with a sample of ~ 2 • 10 4 reconstructed 
decays, which is not far from the capabilities of the present experiments in terms of kaon flux. 
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3 Weak static interaction as the origin of enhancement 



We listed the set of experimental arguments in favour of that a relation between the form 
factors and radiative K decay amplitude takes place. In the following part of the paper we 
would like to show that this relation is not occasional from theoretical point of view. 

In any case, we are trying to address the questions: What are contributions of other loop 
diagrams? What is the origin of the enhancement coefficient g$ in the amplitude jTj)? What is 
origin of the coincidence of the resonance parameters with the kaon decay ones? 

We show below that a reply to all these questions can be the weak static interaction as a 
consequence of Dirac like radiation variables in SM. 



3.1 The Radiation Variables in Standard Model 



As it was shown by Dirac in QED [28], the static interactions in gauge theories are an inevitable 
consequence of the general principles of QFT, including the vacuum postulate. In order to 
obtain a physical vacuum as a state with minimal energy, Dirac eliminated all zero momentum 
fields (with their possible negative contributions into the energy of the system) by solving the 
Gauss constraint and dressing charged fields by the phase factors (see also [29], [30, [31]). This 
elimination leads to the radiation variables and the static interactions in both QED and the 
massive vector boson theory [T4] . 

l(L)/ 



In particular, in QED the radiation gauge-invariant variables A^\A) = 



have propagators J^D^ u (q)J 1 
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J o J o { f 6 <m 

q 2 \ 3 q 2 



[T3] have propagators J^D^ u (q)J u 
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In order to demonstrate the inequivalence between the radiation variables and the Lorentz 



ones, let us consider the electron-positron scattering amplitude T 
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(e + ,e IS^e^e ). One 



can see that the Feynman rules in the radiation gauge give the amplitude in terms of the current 



J 2 
q 2 



ik 
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q 2 



Ji.Jk 



-J 2 



q 2 + is q 2 + ie 



+ 



{qoJo) 2 - (q- j) 
q 2 [q 2 + ie] 



(22) 



This amplitude coincides with the Lorentz gauge one 

1 



T 



q 2 + ie 



J 



(g J - q- J) 

q 2 + ie 



(23) 



when the box terms in Eq. (1221 1 can be eliminated. Thus, the Faddeev equivalence theorem 
[30l is valid, if the currents are conserved 



g Jo - q • J = gJ = 0, 



(24) 



and the box terms are eliminated. It just the case, when the R variables are equivalent to the 
L ones [30] . However, if elementary particles are off their mass-shell (in particular, in bound 
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Figure 4: Axial (a) and vector (b) current contribution into K + — > n + transition 



states) the currents are not conserved^. 

Radiation variables have vacuum as a state with the minimal energy, whereas the Lorentz 
ones lose the vacuum postulate as the time component give the negative contribution into the 
energy. Therefore, Schwinger in [29] ... rejected all Lorentz gauge formulations as unsuited to 
the role of providing the fundamental operator quantization . . . 

Let we believe Schwinger, and consider the massive vector Lagrangian 

£ = -^(W - d^w+^w- - d v W~) + m^w+w- + 

+ [j;w+ + j;w-] 



2V2sm6 w 



in terms of radiation variables p3| W± R = W± + <9 M [1/(M^ - A)]d k W~^ 
In this case, instead of the standard propagator [13] 



-J, 



" q 2 - M? 



w 



(25) 



we have the radiation one 114 



r r 



? + M 2 ) q 



M 2 W 



(26) 



The R-propagator is regular in the limit M\y — > and is well behaved for large momenta. In 
the following we compare two propagators D^ u and D R 



3.2 Weak static interaction as the origin of enhancement 

Let us consider the K + — > n + transition amplitude 

(tt+I -i J dx i dy 4 nx)D^(x~y)r(y)\K + )=z(2n)H\k-p)G EW J:(k 2 ), (27) 

in the first order of the EW perturbation theory in the Fermi coupling constant (El) comparing 
two different W-boson field propagators, the accepted Lorentz (L) propagator (|25j) and the 

x The change of variables R — > L means a change of physical sources. In this case, the off mass-shell L variable 
propagators lose the Coulomb pole forming the Coulomb atoms. The loss of the pole does not mean violation 
of the gauge invariance, because both the variables (R and L) can be defined as the gauge-invariant functionals 
of the initial gauge fields. 
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radiation (R) propagator (1261) .These propagators give the expressions corresponding to the 
diagrams in Fig. [H 



E(^) - E*(fc 2 ) = 2i^ 2 + 2i / ^ f u ' ^ — -, (28) 

J (2tt) 4 (-|g] 2 -M^)[(/c + g) 2 -m2 +ze] 

1 ; 1 ; w J (2vr) 4 (fc + g) 2 - m 2 + *e v ; 

The versions R and L coincide in the case of the axial contribution corresponding to the first 
diagram in Fig. [H and they both reduce to the static interaction contribution because 



1V1 W 



However, in the case of the vector contribution corresponding to the second diagram in Fig. [4] 
the radiation version differs from the Lorentz gauge version (l25ll 2 1. 

In contrast to the Lorentz gauge version fl25l) . two radiation variable diagrams in Fig. S]in 
the rest kaon frame = (k Q , 0, 0, 0) are reduced to the static interaction contribution 

i(2n) 4 S\k-p)G EW E R (k 2 ) = (tt+I -i J dx 4 J (x) A ^ M2 Ux)\K + ) (30) 
with the normal ordering of the pion fields which are at their mass- shelfl so that 



E R (k 2 ) = 2k 2 F 2 



1 + 



M 2 , r dH 



F 2 {2vfJ 2E 7r (l)M 2 v + P 



= 2k 2 F 2 g 8 . (31) 



Here E w (l) = y m 2 + I 2 is the energy of 7r-meson and g$ is the parameter of the enhancement of 
the probability of the axial K + — > n + transition. The pion mass-shell justifies the application 
of the low-energy ChPT [T], where the summation of the chiral series can be considered here 

as the meson form factors OQ El EB1 [ -> [ fMEzWl^iziS^l , 

Using the covariant perturbation theory [32] developed as the series Ju(7 © £) — jZ«(£) + 
F 2 d^ k + J l fijkJ.l,(£) + 0(l 2 ) with respect to quantum fields 7 added to £ as the product 
e *7 e «£ = e «(7®5) ; one can see that the normal ordering 



< 0|7*(a;)7 < '(y)|0 >= 6 ll> N(z), N{z) 



(2^2E~ff) 



where z = x — y, in the product of the currents JZ(j®£) leads to an effective Lagrangian with 
the rule AT = 1/2 



-M w \z\ 

r2 / " - e 



/ d^ 8 (|*|)__[^(£( z ))^(^ +/i.c], 



2 The Faddeev equivalence theorem [30] is not valid, because the vector current J p = Kd^ir — nd^K becomes 
the vertex T M = Kd^D^ — D^d^K , where one of fields is replaced by its propagator OD v — S(x), and d^T^ 7^ 0. 

3 The second integral in ([28| with the term (fco + <7o) 2 really does not depend on k 2 , and it can be removed 
by the mass rotation. 
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where gs(M) = [1 + Yl c 1 N^z)] is series over the multipaticle intermediate states (this sum is 

known as the Volkov superpropagator P||33]). In the limit Mw —* oo, in the lowest order with 
respect to Mw, the dependence of (fed-?!) and the currents on z disappears in the integral of 
the type of 



3 js(\z\)e- M ^ 



d 6 z 



A-n\z\ 



drre r g s (r/M w ) ~ g 8 (0). 



In the next order, the amplitudes K°(K°) — > 7r° arise. Finally, we get the effective Lagrangians 



G, 



£ (A T=i) = ^(0)cose o sine o [(4+i4)(4-<4)-(4+^4)(4-<4)+/i.cj ) (32) 

' " ; ■■" ' " (33) 



£(AT=§) 



— ^ cos 9 C sin 9 C 
v2 



(4 + ^4) (4 



/i.e. 



This result shows that the enhancement can be explained by static vector interaction that 
increases the K + — > n + transition by a factor of g 8 = (fe(0), and yields a new term describing 
the K° — > 7T° transition proportional to g 8 — 1. 

This Lagrangian with the fit parameter g$ = 5 (i.e. g% sin 6c cos — 1) describes the 
nonleptonic decays in satisfactory agreement with experimental data [21 [TT]. 



Conclusions 

We have investigated the low-energy relations between the values of resonance masses in the 
meson form factors, and the differential radiative kaon decay K + — > 7r + e + e~ (/i + /i _ ) rates 
following from the ChPT [UEl [4]. We give non-trivial predictions of muon/electron ratio and 
the effective form factor slope, which are in agreement with the experimental data. 

The high sensitivity of these relations, the low energy status of ChPT, where they arise, 
and the universality of the enhancement coupling constant g% for all kaon-pion weak transition 
amplitudes with the rule of selection AT = 1/2 can be explained by a weak static interaction 
of massive vector bosons in the Hamiltonian approach to SM [T4] . 

The instantaneous character of weak static interaction in the Hamiltonian SM excludes all 
retarded diagram contribution in the effective Chiral Perturbation Theory [8] that destruct the 
form factor structure of the kaon radiative decay rates. The enhancement of kaon-pion tran- 
sition can be considered as a consequence of normal ordering of all pions in the instantaneous 
loop on their mass-shells p 2 = m 2 . The static interaction mechanism of the enhancement of the 
AT = 1/2 transitions predicts the coincidence of the meson form factor resonance parameters 
with the parameters of the radiation kaon decay rates K + — > 7r + e + e~ in satisfactory 
agreement with the experimental data [22l [23]. 

Therefore, the off-mass-shell kaon-pion transition in the radiation weak kaon decays can be 
a good probe of the weak static interactions. 
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A Appendix A: Calculation of K + —> n + l + l~ Decay Width 
A.l The Matrix Element 

The matrix element for the process in Fig. [5] can be obtained by Feynman rules: 

iM{K+ - ir+l+r) = u s (q-)(-ie^)v s '(q + ) i 9 ^- < n + (p)\J em »\K + (k) >, (34) 

after inserting parametrization 

iM(K + -> ir+l+r) = u s {q^){-ieY)v s '{q+) i 9 ^ eF{kp){k + p) u . (35) 

q 2 

To obtain the square root of matrix element one has to sum over spins 
J2\M\ 2 = e4 ' F y - TrUi^- mh/Mp+mp + ky 

s,s' ^ 

Ae 4 \F(kp)\ 2 

= — A [?+a»?-" + q+vV-P- ~ Q+9-9iw ~ m 2 g^ v ](p + fc)^(p + kf 

Ae A \F(kv)\ 2 

= y l 2( l+(P + k).q-(p + k)- q + q-(p + k) 2 - m 2 (p + k) 2 }. (36) 

Next if we define 



q 2 2q+q~ + 2m 2 (p + q-) 2 m 2 + 2pq^ + mf /m^x 2 f mi 



2 



Z M 2 K M£ ' X M 2 K M 2 K ' R \M K ) ' n \M K J ' 

(37) 

where I = e, fi one obtains 

2g4^2 _ . . . . . . 

^|A^| 2 (2,x) = — — \2x + z-2-2n)[-2x- z + 2R + 2n) + z[z-2-2R^ 
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where we used momentum conservation law and relations: 

k = p + q- + q+ = p + q, 

q-q + = %z-2r), 

(k + p) 2 = M 2 K (2 + 2R-z), 
(k-p) 2 = q 2 = M 2 K z, 



q- (p + k) 



M 2 K 



(2x + z-2R-2r), 



M 2 K , 



q+ ( p + k) = ^(2 + 2r-z-2x), 
q(p + k) = M 2 K (l-R). 



A. 2 The decay rate 

The phase volume (in the frame of k — 0) is 



d<5> = 



d 3 q + d 3 q-d 3 p 4 
2e+2e„2e 5 & ~ P ~ ~ ^ 



d A q- = d A q, 



k7+l ' ~ JI ' ( —d\,d>,„rU- - p - q)5((k - q) 2 - ml)5((q - q + f - mj). 



2e 

\q + \de + dQ + 1 



2.2M 



— - o 1 o\ dq 2 A<w5 (q 2 -2q + q); dtt + = 2nd cos 6 + , 
k 2 V / 

4ir2irMl; q + q = q e + - \q\\q+\ cos9 + , 



^4ndq 2 hn = ^ Z —-^ 
8M K H 2 16 



and the decay width is 



4z?2 r 



dT = 



x 



1 2e 4 F 
2M K z 2 

1 dxdz 
(2tt) 5 4 



(2x + z - 2 - 2rj) ( - 2x - z + 2i? + 2rj) + z(z - 2 - 2^ 



7T 2 M|-. 



dT a 2 F 2 



M K dxdz 8nz 2 
After integration the decay width is 



2.r + ; - 2 - 2rj) ( - 2.r - : - 27? + 2r, ) + :( .:. - 2 - 2R 



:(z-2-2R) 



dT a 2 



4r 



where 



A(l, z, R) = (z 2 + 1 + R 2 - 2z - 2Rz - 2R). 
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